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We analyze the generation rates and preparation fidelities of photon triplet states in pulsed cas¬ 
caded parametric down-conversion (PDC) under realistic experimental circumstances. As a model 
system, we assume a monolithically integrated device with negligible interface loss between the 
two consecutive PDC stages. We model the secondary down-conversion process in terms of a lossy 
channel and provide a detailed analysis of noise contributions. Taking variable pump powers into 
account, we estimate the impact of higher-order photon contributions and conversion processes on 
the achievable coincidence probabilities. At mean photon numbers of (m) ~ 0.25 photons per pulse 
behind the first PDC stage, we expect around 4.0 genuine photon triplets per hour. Additionally, 
we discuss fundamental limitations of our model system as well as feasible improvements to the 
detectable photon triplet rate. 

PACS numbers: 42.50.Ar,42.65.Lm,42.65.Wi,42.82.Bq 


I. INTRODUCTION 

Multipartite entanglement has been studied mi and 
experimentally demonstrated (see, for example m in re¬ 
cent years. It has been shown that Greenberger-Horne- 
Zeilinger-states (GHZ) can be used in order to prove the 
nonlocality of quantum physics [4j and, thus, the lack of 
hidden variables m- 

In 2010, Hubei et al. [5: demonstrated the direct 
generation of photonic triplets, which are the funda¬ 
mental states for multi-partite entangled state genera¬ 
tion, using cascaded parametric down-conversion (PDC) 
sources. The same group also verified tripartite entan¬ 
glement in the energy-time regime [9] as well as in the 
polarization degree of freedom m- Their experimental 
results have proven the feasibility of generating heralded 
Bell states with high visibility and fidelity, but with¬ 
out the need of post-selection. Likewise, the violation 
of Clauser-Horne-Shimony-Holt's (CHSH) [TT] inequality 
by over three standard deviations and the Merrnin m 
nu inequality by more than ten standard deviations were 
demonstrated. With this type of cascaded PDC sources, 
hyper-entanglement of photon triplets or heralded GHZ 
states become achievable perspectives. 

In these initial experiments continuous wave lasers 
were employed, while pulsed pump light would provide 
timing information beneficial for source synchronization 
and simplifying the analysis of higher photon number 
contributions. 

Here, we study the photon triplet generation perfor¬ 
mance for the pulsed regime. In detail, we make gen¬ 
eral assumptions concerning the design parameters and 
measurement conditions, and we calculate not only the 
achievable rates of three-fold coincidences, but also show 
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the significant influence of higher-order photon contribu¬ 
tions on the measurement outcome. Furthermore, the 
analysis of the signal-to-noise behavior with real-world 
optical components and detectors is provided. Finally, we 
discuss optimization possibilities for the photon triplet 
source as well as for the measurement apparatus. We find 
that the measurable photon triplet rate can be improved 
applying the latest developments in detector technology 
(see, for example [T5Ul8j l as well as pump lasers with 
high repetition rates. 

The paper is structured as follows: In Section [Hj we 
briefly describe the considered model system, which is 
followed by the estimation of primary PDC photon num¬ 
ber output in Section HI Sections IV and |V| describe the 
secondary PDC process in terms of a lossy channel and 
the spectral splitting of the generated secondary photon 
pairs, respectively. In Sect ion [VT| we analyze the two- and 
three-fold coincidence probabilities, and we conclude in 
Section IVIII 


II. DESCRIPTION OF THE MODEL SYSTEM 

Our theoretical and experimental considerations are 
inspired by the idea to generate photon triplets with cas¬ 
caded PDC processes, as shown by Hubei et al. 0. As 
a model system we consider the monolithic guided-wave 
device depicted in Fig. [l] The input waveguide, assumed 
to be fabricated by titanium-diffusion for low loss HU, is 
periodically poled with two different grating periods in 
consecutive areas. In the first poled region, nondegen¬ 
erate (primary) photon pairs at around A s i ~ 790.5 nrn 
(signal 1) and An ~ 1625 nm (idler 1) are generated in a 
type-0 parametric down-conversion (PDC) process from 
a pulsed pump at 532 nm. 

In order to separate the primary pair of photons in a 
spatio-spectral manner, we consider a passive wavelength 
division multiplexer (WDM) between the two poled re- 
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gions. Idler 1 photons at the longer wavelength are trans¬ 
ferred to an adjacent waveguide by evanescent field cou¬ 
pling with efficiency Coupler 41 , and they pass on-chip op¬ 
tical losses with efficiency r/\y g.u- The WDM is wave¬ 
length selective in the telecom region for transverse mag¬ 
netic (TM) polarization. 

The short-wavelength signal photons, generated in the 
fundamental spatial mode, remain in the pump chan¬ 
nel with efficiency //coupler,si due to strong mode confine¬ 
ment, and they pass optical loss along the WDM and the 
straight waveguide with efficiency p wc,si- 

The aforementioned components of the model system 
have been experimentally tested in {20 GO] , and we as¬ 
sume similar linear and nonlinear optical properties here. 

In the secondary periodic poling, primary signal pho¬ 
tons are predominantly present and decay with the nom¬ 
inal conversion probability Pp£,“ 2 into secondary photon 
pairs (signal 2 and idler 2 ), which are distributed around 
the degeneracy wavelength of A Sj 2 = Ai ,2 = 1581 nm. 
This has to be in accordance to the rules of energy and 
momentum conservation (phase-matching) for the sec¬ 
ondary type-0 PDC process. We assume that detuning 
of the secondary down-conversion source allows for non¬ 
degenerate operation such that the overall cascaded de¬ 
cay of a green pump photon into three telecom photons 
is described by 

532 nm —> 1625 nm + 1551 nm + 1611 nm (1) 
pump —>• idler 1 + signal 2 + idler 2 . 

In the case of primary signal photons we expect loss 
figures of pwG.si < 0.3dB/cm and ?/wG,ii = Pwg,s2 = 
?/wG,i 2 < 0.08dB/cm for telecom photons. All generated 
telecom photons can exit the device through an end-face 
coating with efficiency pef > 0.99. 

Although our model system is an integrated device, 
the following considerations apply also to bulk and hybrid 
cascaded parametric down-conversion sources, if they are 
pumped with a pulsed laser system. Figure [T] shows the 
principle measurement scheme as well as the variables 
used in this work. 


III. ESTIMATION OF THE MEAN PHOTON 

NUMBER IN THE PRIMARY PDC PROCESS 

In this section we model the outcome of the primary 
parametric down-conversion process in terms of the click 
probability for idler 1 photons at binary detectors. We 
assume the signal and idler photon number per optical 
pulse to be Poisson-distributed, since type-0 PDC pro¬ 
cesses are typically multi-mode in the spectral domain 

EB- 

The density vector element, which describes the prob¬ 
ability that an m-photon state contributes to the average 
photon number (m) per signal or idler pulse, is given by 


Depending on the pump power, i. e., the photon num¬ 
ber density per pump pulse, we can deduce the mean 
photon number per output pulse fairly well by measur¬ 
ing idler 1 clicks events with a binary free-running detec¬ 
tor. This is only valid under the precondition that the 
detector does not saturate due to dead-time effects. 

For the (lossless) case of an arbitrary m-photon input 
state impinging the idler 1 detector, we can model the 
click response of the detector. If no photons arrive (m = 
0), only noise counts can influence the zero-photon click 
probability, which reads 

Pi 1 (TO = 0) = Pnoise.il ■ (3) 

The noise count probability of a free-running detector, 
Pnoise.ii j can be deduced from the specified dark count 
rate Rdc,n, and from the count rate of additional black- 
body photons, -Rbb- The latter becomes significant, 
if we heat up the cascaded PDC source to operation 
points around 170°C. We relate the sum of the two rates, 
-Rnoise.ii; to appropriate timing units in order to get the 
noise probability Pnoise.ii- For pulsed systems, the pulse 
duration, in combination with the joint timing jitter of 
the involved devices, appears to be a reasonable timing 
choice. 

Second, when we know that there is exactly one photon 
impinging to the detector (m = 1), a click can occur 
(logical OR) due to a noise count or due to the primary 
idler photon itself with efficiency 77 Det.ii or both, i. e., 
noise counts ’blind’ the detector for the incoming photon. 
Therefore, the are treated as stochastically independent, 
and we get the click probability for the one-photon state: 

P[[ (m = 1) = Pnoise.il T PDet.il Pnoise.il * PDet.il- (4) 

As soon as there are (higher-order) m-photon states ar¬ 
riving at the detector, the latter will give only one click, 
since it is not photon-number resolving. The correspond¬ 
ing click probability is deduced as the counter event of 
the case, where the detector neither clicks due to noise 
events (that is p = 1 — p no ise.ii); nor due to any of the 
impinging photons (i. e., p' = (1 — pDet.ii)" 1 )- This corre¬ 
sponds to a logical NAND operation to the negated zero- 
and one-photon click probabilities, and we write 

P'il ( m ) = 1 - (1 - Pnoise.il) ' (1 ^ PDet.il)” 1 - (5) 

Since this formula is consistent with equations [3] and [4] for 
the cases m = 0 and m = 1, respectively, we generalize 
equation [5] for arbitrary m-photon states with m > 0. 

We expect primary idler photons to undergo optical 
loss of 1 — pH on their path from the point of generation 
to the point of being detected. In detail, the transmit¬ 
tance pii can be calculated as the product of individual 
transmittances and efficiencies of lossy optical elements 
and the detector (logical AND), that is 

Pil = Pint.l • PEF • POC.il ' PDet.il- (6) 

Herein pint.i = Pcoupier.il ■ PwG.ii denotes the intrinsic 
transmittance of our model device for idler 1 photons 
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FIG. 1. Monolithic source and measurement schematic for cascaded parametric down-conversion. In the periodically poled 
area I, pulsed pump photons decay to nondegenerate photon pairs, si and il, sharing the same mean photon number per 
pulse ( m). The primary photon pair splits up spatio-spectrally at the integrated wavelength division multiplexer (WDM). 
Short-wavelength si photons remain in the original waveguide and decay to time-correlated secondary photon pairs, s2 and 
i2 (with common mean photon number per pulse (k) <C (m )), in the periodically poled area II. Effectively, we convert one 
pump photon into three telecom photons, which are directed to individual detectors and undergo coincidence measurements, 
as described in detail in the text. 


behind the first PDC stage. The factor 770 c,il is the 
transmittance of supplementary optical components in 
the idler 1 measurement arm, and the detector efficiency 
is still labeled ?7Det,ii- With this in mind, we substitute 
equation [5] by the loss-dependent m-photon click proba¬ 
bility 

Pi (m) = 1 - (1 - _Pnoise.il )(1 - Vil) m - (7) 

Third, we calculate the overall click probability 
Pi 2 ((to)) for an ensemble measurement of optical pulses, 
carrying Poisson-distributed photons with mean photon 
number (to). We weight the individual m-photon con¬ 
tributions with their occupation probabilities p rn and get 
for the overall click probability 


one-photon state. Furthermore, there is a linear relation 
between Pi ((to)) and (m) for (m) <C 1, where only zero- 
and one-photon components play a predominant role. 

Next, we model the response of a binary detector for 
a real-world scenario. If p\\ is properly known, we can 
calculate the click rate of a free-running detector in the 
idler arm to be approximately 

Pil ~ Pnoise,il T Prep ' Pil ((to)) , ( 9 ) 

where P rep is the pump laser repetition rate. 

Comparing the measured click rates with the calcu¬ 
lated ones, we can deduce the mean photon number of 
the first PDC process by solving equation [9] for Pi ((m)) 
and get 


OO 

Pil ((TO)) = y~) [1 - (1 - Pnoise,il)(l - Vil) m ] ' Pm- (8) 
m> 0 

We find that Pi ((to)) is equal to the noise count prob¬ 
ability p n oise,ii for (to) —> 0. Contrarily, for large mean 
photon numbers (to) 1, Pi ((to)) asymptotically ap¬ 
proaches unit click probability. Another consequence of 
equation [8] is its linearity with respect to (to) for low 
arm efficiencies r/n -C 1. This means that an m-photon 
state will yield a click about m-times more often than a 


Pi ((m)) = 


Pil — R 


noise.il 


P 


rep 


Pi: 

p" 


( 10 ) 


rep 


For small arm efficiencies <C 1, we can infer that 
the mean photon number per pulse is approximately 


( to ) 


Pil 


( 11 ) 


Vil ' Prep 

Note that we assumed in the measured click rates 
in equation 10 to be at moderate levels, where neither 
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dead-time effects of the detector nor its saturation 
become significant. The validity of the latter condition 
depends on the efficiency pn. However, the overall count 
rate should exceed the noise count rate significantly, i. 
e., R\\ -ffnoise.iij for the approximations in equations 
flOl and flOl to hold. 


IV. SECONDARY PDC PROCESS 

Accessing the mean photon number of the first PDC 
process by measuring primary idler clicks also allows us 
to predict the outcome of the secondary PDC process, 
since the primary photons share temporal correlations 
and have the same mean photon number per pulse at the 
point of their generation. Thus, we are able to estimate 
the contributions of individual m-photon states in an en¬ 
semble of optical pulses. We can be sure that, at the exit 
of our primary PDC stage, the photon number distribu¬ 
tion of the signal photons is not only Poisson-like, but 
also equal to the corresponding idler distribution, if we 
neglect different loss values inside the chip. 

In the following, we model the transfer of signal 1 pho¬ 
tons towards the secondary conversion stage as well as 
the second PDC process itself in terms of a lossy channel 
with subsequent deterministic pair generation. An in¬ 
coming m-photon state behaves as if it is reflected (and 
lost), when no conversion occurs, but gets transmitted 
and duplicated, when a conversion takes place. After 
conversion, we treat the generated pair photons as indi¬ 
vidual but time-correlated entities. 

We also assume that the spectral characteristics of sec¬ 
ondary photon pairs can be tuned such that determinis¬ 
tic splitting is achieved. Experimentally, this could be 
realized by implementing a long-pass (LP) filter and by 
driving the secondary PDC at nondegenerate photon pair 
emission. 

The linear optical device efficiency, 

^?Dev,sl ~ ^coupler,si ' ^WG,sl> (12) 

describes the transmission of signal 1 photons through 
the integrated WDM towards the secondary PDC stage, 
while ?7coupier,si is the efficiency that the signal photon 
remains in the original coupler arm. The term ?7wG,si 
denotes the waveguide transmittance for the signal 1 pho¬ 
ton from the average point of its generation to the average 
point of its decay to a secondary photon pair. 

The secondary PDC process takes place with the nom¬ 
inal conversion efficiency per signal 1 photon Pp£,“ 2) 
which includes the spectral overlap of the primary sig¬ 
nal photon mode with the phase-matched pump mode 
of the secondary process. Note that the mutual spatio- 
spectral mode compatibility between signal 1 photons 
and the pump input of the second PDC process is impor¬ 
tant, but we restrict ourselves to include this precondi¬ 
tion to PpD™ 2 - We multiply the probabilities/efficiencies 


accordingly, since they describe stochastically indepen¬ 
dent events, and we get for the internal conversion prob¬ 
ability per incoming signal 1 photon in the second PDC 
process 

PpDC,2 = ?7Dev,sl • PpDC,2i (13) 

which acts similar to a “loss factor” on primary signal 
photons. Note that ?7Dev,si can be replaced by the trans¬ 
mittance of any interface as long as the considered model 
system is a bulk or hybrid cascaded parametric down- 
conversion source. 

We write the resulting photon number transformation 
in a matrix representation: 

Pk = LpDC,2 • Pm- (14) 

The vector p m comprises the incoming photon number 
occupation probabilities for different m-photon states ac¬ 
cording to equation [2] The elements 

Lkm = ^ ^ ^PDC,2 (1 — PpDC,2) m k , m>k (15) 

contribute to the loss matrix Lpdc ,2 and describe the 
probability that k secondary photon pairs are generated, 
given an m-photon state in the signal 1 (and idler 1 ) 
mode is present behind the first PDC stage. 

This leads to individual photon number contributions 
Pk of the time-correlated signal 2 and idler 2 photons. 
We write 

Pk = ^2 ( j^J -PpDC,2 (! “ -PpDC,2) m k • Pm- (16) 

In the limit of unit internal conversion efficiency 
PpDC ,2 = 1, we expect pk = p m to hold true, whereas 
real-world conditions lead to drastically decreased mean 
photon numbers per pulse for secondary photons, i. e., 
(k) -C (m), since nominal conversion efficiencies are 
typically of the order of 10 -5 < Pppc y2 — 10 -10 pairs 
per input photon. 


V. SPECTRAL SPLITTING OF SECONDARY 
PHOTONS 

In order to perform coincidence measurements be¬ 
tween secondary photons, we have to split them spatio- 
spectrally. A realistic case for quasi-deterministic sepa¬ 
ration of nondegenerate secondary photon pairs is given, 
if we use, for example, a long-pass filter with a steep cut- 
on edge as the splitting element. Additional band-pass 
filters could be inserted in both, the reflected and the 
transmitted beam path. This provides the reduction of 
noise events, which are related to the blackbody emission 
of a heated cascaded parametric down-conversion source. 

We combine the wavelength dependent transmission at 
supplementary optical elements (filters, fiber-couplings et 
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cetera) with the splitting behavior of the long-pass filter 
and with the detection efficiency in the two individual 
measurement arms: 

Vs2 (As2) = -RlP (As2) • ? ?OC,s2 (A s 2) ' VDet,s2 (A s 2) • ?7EF (A s2 ) 

(17) 

and 

Vi2 (Ai 2 ) = Tlp (Ai 2 ) • ??OC,i2 (A12) • Wet,i2 (Atf) ’ VEF (Ai 2 ) • 

(18) 

Herein, the optical components provide transmittances 
770 c (A) in the respective beam path. Note that tjef (A; 2 ) 
describes the wavelength-dependent transmittance of 
telecom photons at the device’s end-facet. Individual 
detector efficiencies are labeled with r/Det (A), the wave¬ 
length dependence of which must be taken seriously, if 


common InGaAs-based detectors are used. 


The aforementioned considerations allow us to antic¬ 
ipate the click probabilities of binary detectors in the 
respective measurement arm. The derivation is similar 
to equation [8j but we take the significant changes of the 
photon number distribution, caused by the second PDC 
process (see equation 16), into account. Additionally, we 
pay attention to the spectral dependence of the transmit¬ 
tances in equations EH and [18] 


A click event in one of the free-running detectors is 
given as the counter event of having neither a dark count, 
nor signal/idler photons, respectively. For the click prob¬ 
abilities of single detection events in the respective sec¬ 
ondary arm we write 


OO 

P&2/\2 {( m )) — _ (l — Pnoise,s2/i2) [l — Vs2/i2 (^s2/i2)] } ' Pk 

k> 0 

00 00 / \ 

= ^ - Pnoise,s2/i2) [l ~ Vs2/i2 (A s2 /i2)] ^ ^ x ( jfc ) ^PDC,2 (1 ~ ftDC,2) m * Pm- 

k> 0 ' 


rri>k 


(19) 

( 20 ) 


We consider p n oise,s 2 /i 2 to be the individual noise count 
probabilities of the signal 2 and idler 2 detectors. With 
the realistic assumption of Pp£,™ 2 ~ -Ppdc ,2 ~ 10 -7 , 
we can conclude that all (m > l)-photon states of the 
primary PDC process will be converted mainly to k = 0 
and few k = 1 photon contributions in the secondary 
PDC process. 

In the following, we assume that the wavelength de¬ 
pendent arm efficiencies can be described by analyti¬ 
cal expressions, being continuous for reasonable inter¬ 
vals around the expected signal/idler wavelengths, i. e., 
AW" 2 < A s2 /i 2 < A“^. Proper intervals given by the 
filter bandwidths of, say, fiber-optic band-pass filters. 
Thus, we write for the overall arm efficiencies 77 tot , which 
are experimentally accessible: 


•jg* 


f Vs2 (A S 2) dA 


\ max \ min 

A s2 A s2 


and 


( 21 ) 




/ r] i2 (A i2 ) dA 


\ max \ min 

A i2 A i2 


( 22 ) 


the click probabilities for secondary PDC photons: 

OO 

-fs2((^)) ~ Pnoise,s2 “I” Ps2 ^ ^ ^ ‘ -^PDC,2 * Pm (23) 

ra> 1 

i „^tot ngGIl /O/lA 

— Pnoise,s2 + 7 7 s 2 ' * \ PDC,2 (24) 

— Pnoise,s2 T ?7 S 2 ’ PpE>C,2 * (^) • (2b) 

and 

OO 

P’i2 ((^)) ~ Pnoise,i2 “b Tj\ 2 ^ ) TTl * PpDC,2 * Pm (26) 

m> 1 

= Pnoise,i2 + Vtt* ' ^PDC,2 (27) 

= Pnoise,i2 + ’ -PpDC,2 ' (m). (28) 

In equations [23] to [28] we notice, that primarily gener¬ 
ated higher-order ?n-photon states will each contribute to 
the secondary conversion process approximately m-times 
more often than one-photon states. Note that the term 
TpDC 2 labels the effective pair generation probability per 
optical pulse, whereas Ppdc ,2 is the pair generation prob¬ 
ability per incoming pump photon. 

We infer that individual click probabilities depend lin¬ 
early not only on p no i S e,i 2 , -Ppdc, 2 , and on r 7s2/i2 ’ but 
also on the mean photon number of the primary PDC 
output,(m), since 


We furthermore approximate (1 — PpDC, 2 ) m k ~ 1 for 
reasonable internal conversion efficiencies. Additionally, 
we expect p no i se , S 2 < vl? and p no ise,i 2 < . 

With all these assumptions and preconditions, the 
evaluation of equation 20 for both photon species yields 


^2 rn ■ p m = (to) . (29) 

ra> 1 

This means that four parameters have an influence on 
the single-click probabilities P((m)). While the noise 








6 


probability has to be kept as low as possible for good 
signal-to-noise ratios (SNR), the efficiency 77 * 2/12 should 
be as high as possible. 

Increasing Ppdc ,2 to values higher than 10 -6 is tech¬ 
nologically hard in lithium niobate waveguide structures. 
This fact implies, for experiments with commercially 
available equipment, that the noise count probability is 
typically of the same order as Ppdc, 2 , and we will not be 
able to identify photon triplets with good signal-to-noise 
ratios just by measuring the generation of secondarily 
generated signal or idler photons. Instead, we will make 
use of the temporal correlations of the triplet photons 
and perform coincidence measurements, as described in 
the following. This will lift the genuine triplets above the 
inevitable accidental background contributions. 

Note that equations [23] and [26] reflect what we would 
also expect for direct pumping of the secondary PDC pro¬ 
cess with attenuated laser pulses instead of single photons 
from the primary PDC stage. 


VI. COINCIDENCE MEASUREMENTS 

For the evaluation of two- and three-fold coincidences, 
we look at the possible outcome options: in practice, we 
must accept that clicks at free-running binary detectors 
will stem from noise contributions OR genuine PDC pho¬ 
tons. 

We distinguish between eight cases as shown in Fig. [2] 
where only case A represents a genuine photon triplet. 
All other options must not be discarded, but have to be 
taken into account as uncorrelated noise contributions. 


A. Two-fold coincidences of secondary PDC 
photons 


The conditional detection probabilities (or Klyshko effi¬ 
ciencies 77 K [22]) are seemingly independent of the single 
photon pair generation probability, and they represent 
the respective overall path efficiencies. 

It has been shown experimentally that the coincidence 
probability scales super-linear for pump pulses, which 
contain several trillions of photons, while the single count 
probabilities scale almost linear with increasing pump 
powers m ■ This effect seems to be contradictory to 
the aforementioned considerations, but so far we implied 
the generation of twin photons only, while we neglected 
higher-order photon contributions appearing at higher 
pump powers. 

In our case, we recognize the similarities of equations 
[24] and [27] and equations [30] and [3l] respectively. Hence, 
we infer the probability for the coincidental detection of 


the time-correlated secondary photons from equation 32 
and we write for the noiseless detection: 


P 2 c ”fo ld (s2 D i2, (m)) = ■ P|“ c>2 (35) 

OO 

= vlTvPf 1 m ' P pdc ,2 • M.36) 

m> 1 

= rfsfrfiT ■ -Ppdc ,2 • (m). (37) 

This means that the coincidence probability of secondary 
pairs scales linearly with the mean photon number of the 
primary process. 

We estimate the noise-related accidental coincidences 
to be 


A-fold ( s2 n i2 > ( m )) = vl? ■ fl?DC,2 • (m) • Pnoise,i2 (38) 

' -^PDC,2 ’ {jTl) ‘ Pnoise,s2 
TPnoise,s2 ' Pnoise,i2; 

and calculate the signal-to-noise ratio for secondary co¬ 
incidence detection: 


Parametric down-conversion processes in the pulsed 
regime generate time-correlated pairs of photons, which 
must not be treated as independent entities. In principle, 
the individual and the coincidental detection probabili¬ 
ties both scale linearly with the pair generation probabil¬ 
ity Ppdc- We write for the signal and idler probabilities 

P s = 77 s • Ppdc (30) 

Pi =T7i-PpDC) (31) 

and for the coincidence probability of correlated events 
we have 

P (s (~l i) = P; • P s (s|i) = P s • Pi (i|s) = 77 s ? 7 i • P PD c, (32) 

where r] s and 771 represent the overall efficiencies in the 
respective detection arms. The detection of one photon, 
given that the detection of its twin already occurred, is 
expressed by the conditional probabilities 

p s (s|i) = P ^ ^ = 77 k,s = 77 s (33) 

-*i 

p i (i|s) = P ^ ^ = 77 k, i = 7 ?i (34) 

S 


5'VPpDC,2 


Pf-foid (32 012 , (to)) 
^ 2 n -ioid (s 2 n i 2 , (to)) ' 


(39) 


This result indicates the feasibility to verify photon 
triplet generation by detecting coincidences only of sec¬ 
ondary PDC photons, as long as the value of SNRy>t>c ,2 
is sufficiently high. However, for an imperfect triplet 
source, we have to consider the spatio-spectral split¬ 
ting behind the first PDC stage to be nondeterminis- 
tic. As long as implemented filters for secondary PDC 
wavelengths provide sufficient extinction at primary idler 
wavelengths, we can assume that the probability of de¬ 
tecting idler 1 photons in one or both of the secondary 
measurement arms is negligible and, thus, does not con¬ 
tribute to accidental coincidences of the secondary pho¬ 
ton detection. 

The situation is different and more critical in the case 
of higher-order mode combinations in the first PDC pro¬ 
cess of our model system. These can occur, if the 532 nm 
pump and/or the signal 1 ’ photons are exited in, e. g., the 
TM 0 1 mode and propagating along the waveguide struc¬ 
ture. Due to dispersion and different phase-matching 
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FIG. 2. Possible outcome of two- and three-fold coincidence measurements between the primary and secondary detection 
events. Only case A represents genuine photon triplets preparation, whereas all other options denote noise-related and accidental 
coincidences. The dashed line marks the spatial separation of events with high (top) and low (bottom) detection probabilities. 


conditions for higher-order PDC generation, we must as¬ 
sume that the corresponding idler T photons, although 
in the fundamental TMqq mode, can have wavelengths 
similar to the secondary PDC photons, i. e., Am ~ An, 
which the on-chip coupler might not be optimized for. 
Given this case, the parasitic idler 1’ photons could pass 
the WDM through the original arm (with transmittance 
^coupler,ii' = 1 - ^coupler,ii'(Am). They could also pass 
subsequent optical elements and band-pass filters to¬ 
wards the respective detector. Their contribution to sec¬ 
ondary detection events is given by the probability, that 
at least one of the accidental idler 1’ photon survives in 
the wrong measurement arm: 


we deduce a linear increase of P{^ a2 and P i with the 
pump power, i. e., the average photon number (m / ), for 
small values of CSpier.ii- 

Falsely directed idler 1’ photons have an impact also on 
the coincidence click probability of the secondary process. 
We find that accidental coincidences occur with the joint 
probability 

K Sc 1 ' (s2 n 12, (m» = P^ s2 ■ p$% i2 (42) 

+PiV% S2 • Pi 2 ((m)) 

+PwU i2 • Ps 2 ((m)) 

^ pace pace 
~ Ml'->-s2 ' Ml'— m2> 


and 


Pn'%s2 = V s 2 ’ Vwg,h' ■ Vef 

OO 

x 1 1 - 1 1 - ^ 
m> 0 

P?v%i 2 = vtT ■ ?fwG,ii' • Vef 

OO 

^ [i - (i - V, 


(40) 


ace 
coupler, 


il')"*] -P'r, 


(41) 


X 


m> 0 


ace 
couple: 


r,il' / 


Pm ‘ 


For simplicity we assume that the internal conversion 
efficiency for the primary PDC process is equal for all 
spatial mode combinations. Thus, photon number occu¬ 
pation vector elements p' m = p m , identical to the ones of 
the fundamental PDC process, contribute to the formu¬ 
las above. In practice, the conversion efficiencies strongly 
depend on the waveguide properties (i. e., the effective 
refractive indices of the involved modes) as well as on 
the possibility for quasi-selective excitation of different 
spatial pump modes. 


We notice in equations 40 and 41 that the individual 
click probabilities for accidental idler 1’ photons in one 
of the secondary measurement arms strongly depend not 
only on the properties of the integrated WDM structure, 
but also on the initial pump power and conversion effi¬ 
ciencies for primary higher-order PDC processes, both 
determining the individual photon number occupation 
densities p' m . With the same arguments as in Section III 


which includes coincidences of accidental idler 1’ pho¬ 
tons with noise contributions. For the approximation, 
we accounted for the typically low values of Ppdc ,2 •C 
Cupier.il'- However, if we faithfully assume Cupier.il' < 
1 , we can approximate the coincidences-to accidental- 
ratios (CAR) for the secondary photon detection accord¬ 
ingly: 


Ch4f?pDC,2 


P c c 0 T n r c (s2ni2,(m)) 
CXf (s2 PI i2, (m» 


Ppdc,2 • {m) 


’JwG.ii' ' Vef ' (v 


acc 

coupler, il' 


( to ')) 2 


(43) 


We notice that the coincidences-to-accidentals ratio de¬ 
creases rapidly for higher pump powers, which denotes 
a fundamental limitation of an imperfect integrated de¬ 
vice. Practical ways to overcome the decreasing fidelity 
is, consequently, to provide the suppression by higher- 
order PDC processes in the first stage by proper waveg¬ 
uide engineering. Linearly tapered waveguides, for exam¬ 
ple, can prevent higher-order pump modes at 532 nm to 
be guided. Likewise, tapered structures in between the 
integrated coupler and the secondary poled region can 
provide primary signal to be guided only in the funda¬ 
mental spatial mode, whereas higher-order spatial modes 
at signal 1’ wavelengths as well as telecom photons from 
the first PDC process scatter to the substrate due to the 
shallower refractive index profile of the waveguide behind 
tapering. These technological countermeasures can also 








reduce the necessity for strong and narrow-band filtering 
in the signal 2 and idler 2 measurement arms. 

B. Three-fold coincidences 

Although the coincidental detection of photon pairs 
generated in the second PDC process (provided that val¬ 
ues of SNRpdc ,2 and CARpp,c ,2 are sufficiently high) 
can be seen as the sufficient condition to verify the gen¬ 
eration of photon triplets, the full prove in terms of a 
necessary condition is given only, if we also detect the 


corresponding idler photon from the first PDC process 
in the dedicated measurement arm. 

We rely on a measurement scheme, where all detector 
outcomes are fed into a multi-channel time-tagging unit, 
and we post-select the three-fold coincidences. In this 
configuration, we benefit from having access to acciden¬ 
tal events by analyzing coincidences between neighboring 
pulses |20j . which appear at multiple integers of the in¬ 
verse repetition time of the pulsed pump laser. 

The three-fold coincidence probability of detecting 
time-correlated photons, which include genuine photon 
triplets, is given by 


P (i2 fl s2 D il, (m)) = ^ m • Ppdc ,2 ■ (1 — PpDC, 2 ) m 1 • [1 — (1 — Vii)™} ■ Pm 

m> 1 

oo 

~ vl? ■ vt? ■ Ppdc ,2 ■ ^2 m • [1 - (1 - T]i i) m ] • p m . 

m> 1 


The photon triplet detection probability is given by 

P (i2 n s2 n il, m= 1) « rffi • vlT ■ P PD c ,2 ■ Vn ■ Pi- 
Thus, we can rewrite equation |45| as 


P (i2 fl s2 fl il, (m)) 


P (i2 fl s2 n il, to = 1) 


1 + V. m ■ [1 - (1 - va) m ] ■ ~ 


i>2 


Pi 


(44) 

(45) 


(46) 


(47) 


if we neglect noise contributions and assume, again, 
(1 — PpDC, 2 ) m-fe ~ 1. We notice the influence of 
higher-order photons on the overall three-fold coinci¬ 
dence probability P (i2 fl s2 n il, (m)). Thus, we calcu¬ 
late the factor r((m)) as the ratio of three-fold coinci¬ 
dences, that stem from (m)-photon contributions, and 


(to = l)-photon events using equations 46 and 47 


P (i2 fl s2 fl il, (m)) 

P (i2 n s2 n il, to = 1) 

OO 

m>2 Pl 


(48) 


This indicates that higher-order photons in the first PDC 
stage have a nonlinear impact on the three-fold coinci¬ 
dences. We note that the ratio ^ is linear only for 
to = 2. For to > 3 we expect quasi-linear behavior only 
for low mean photon numbers (to) <C 1 , whereas the 
higher-order photon impact is super-linear otherwise. 

The three-fold coincidence probability increase is de¬ 
picted in Fig. [3j where we plotted r((m)) color-coded 
and dependent on the mean photon number (to) and the 
primary idler arm efficiency pn. We deduce from the 
adjacent graphs that, on the one hand, the idler 1 arm 
efficiency impacts r ((to)) quasi-linearly for (to) ss 0.25. 
On the other hand, r ((to)) increases super-linearly with 


(to), which limits the scalability of a pulsed source fun¬ 
damentally in terms of the injected pump power. Note 
that the inverse of r ((to)) is a measure for the ratio, with 
which genuine photon triplets contribute to the overall 
three-fold coincidence probability. This ratio decreases 
accordingly with increasing pump powers at the cascaded 
PDC source input. 

Consequently, we treat three-fold coincidences, which 
involve higher-order photon contributions of the first 
PDC process, as accidentals. We have experimental ac¬ 
cess to the approximate value of r ((to)) by analyzing 
three-fold coincidences, which include neighboring pulses 
of the idler 1 detection at multiple integers of R rep - The 
probability of detecting three-fold coincidences in two 
consecutive pulses is equal to the probability of register¬ 
ing a three-fold coincidence between a secondary pho¬ 
ton pair and a two-plroton idler 1 pulse. Thus, for 
mean photon numbers (to) < 1 , where the higher-order 
photon states are dominated by two-photon contribu¬ 
tions, we can approximate the accidentals in this manner. 
Note that this method is only possible in pulsed sys¬ 
tems, whereas continuous-wave realizations would have 
to rely either directly on photon-number-resolving detec¬ 
tors, or on measurements of the second-order autocorre¬ 
lation function of primary idler photons. 

We calculate the achievable coincidences-to-accidentals 
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r(<m>) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


mean photon number <m> 


FIG. 3. Increase of the three-fold coincidence click probability dependent on the mean photon number and the idler arm 
efficiency of the primary PDC stage. The adjacent graphs show the individual dependencies for the dotted cross-cut lines, which 
represent reasonable values of rjn and (m). In general, we notice a nonlinear dependence of r((m)) on both axis parameters, 
which represents a fundamental limitation on the generation and verification of genuine photon triplets in a pulsed system in 
conjunction with binary detectors. 


ratio as 


we derive for the expectable photon triplet rate: 


C^4i?3_fold 


P (i2 n s2 fl il, to = 1) 

oo 

P (i2 fl s2 n il, to) 

m>2 

l 


r ((to)) — 1 ’ 


(49) 


which does not include noise-related contributions to the 
three-fold coincidences. Note that the coincidences-to- 
accidentals ratio has values CAR 3 _ fold < 1 for r ((m)) > 
2. This means, in turn, that less than half of the three¬ 
fold coincidences stem from genuine photon triplets. We 
estimate C AR 3 _f 0 \d ss 3.3 for reasonably low pump pow¬ 
ers, i. e., (to) = 0.25, and an idler 1 arm efficiency of 

m = 0 - 117 . 

The question, how many genuine photon triplet states 
we can expect to measure, is of interest now. From 
the approximation in equation [46] and the efficiencies as¬ 
sumed for a realistic experimental setup (see Appendix), 


^triplet = P (i2 ("I s2 n il, to = 1) • i? rep (50) 

= 7 7i2 >t ' vlT ' -PpDC,2 • Vil ■ Pi ■ Rlep 
~ WiS* ' Vl? ’ -FPDC,2 ’ Ril 
= 4.04 h -1 , 

where a mean photon number of (m) = 0.25 per op¬ 
tical pulse in the first PDC process and a reasonable 
pump laser repetition rate of R rep = 10 MHz are in¬ 
cluded. From the result in equation [50] together with 
the individual measurement arm efficiencies, we derive 
a photon triplet generation rate inside the monolithic 
model system of -Rj)7pi e t ~ 1765 h -1 . The theoretical 
generation rate is dependent on the repetition rate of the 
pump laser and the internal conversion efficiency of the 
secondary PDC process. Thus, it can be increased with 
appropriate technological improvements of the influenc¬ 
ing parameters. 

We assume picosecond pump pulses, and we expect the 
three-fold coincidences to be distributed over a temporal 
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width of ~ 1 ns, which includes the joint timing jitter 
of the detection apparatus. By analogy to the case of 
secondary photon pair coincidences, as discussed in Sec¬ 
tion VIA we derive the signal-to-noise ratio SNR 3 _{ a \ d 
by relating the probability of generating three-fold co¬ 
incidences by cascaded PDC, i. e., P (i2 fl s2 n il, (m)), 


to the probability of detecting noise-contributed three¬ 
fold coincidences. A valid approximation of SN i? 3 _f 0 id 
includes three-fold coincidences, where the idler 1 detec¬ 
tion events (PDC photons and noise contributions, see 
equation [8]) as well as the two-fold coincidences of the 
secondary photon detection are involved: 


P 3-fo e id ( i2 n s2 C il, (to)) « P 2 C -foid (s2 n i2 , (to)) • Pnoise.il + Pf-foi id (s2 O i2, (r 






= ■ PpDC.2 • (to) • Pnoise.il + P£-Toid ( s2 0 i2 > ( m )) ' t 1 “ ^ 

m> 0 


(51) 


\ m l 

Pil) I ' Pm 


P 2-foU ( s2 n i2 > ( TO )) • t 1 “ ( 1 “ ' P- 

m> 0 


The approximations in the formula are valid for the 
realistic case, where SNRpnc 2 • Pnoise ii •C Pii holds 
and P^fold (s2 D i2 , (to)) « > 2 c °£ ld (s2 n i2 , (to)) « 
Pii ((to)). For the upper bound of the signal-to-noise 
ratio of three-fold coincidences, we calculate 


SNR 3 _ lold 


P (i2 fl s2 fl il, (to)) 

P 3-fo e id (i2ns2nil,(TO)) 

P (i2 fl s2 fl il, (to)) 
^ 2 -iold (s2 n i2, (to)) • Pi! ((to)) ' 


We see that SNR3 — fold depends strongly on the proba¬ 
bility of noise-related two-fold coincidences and, in turn, 
on the It is important to mention that this measure does 
not represent the signal-to-noise ratio of genuine photon 
triplets with respect to the noise-related background, but 
it rather depends on the overall three-fold coincidence 
probability. Thus, it is highly recommended to verify the 
generation of photon triplets by post-processing analy¬ 
sis of the signal-to-noise ratio and the coincidences-to- 
accidentals ratio in pulsed cascaded parametric down- 
conversion processes. 


VII. CONCLUSION AND OUTLOOK 

We presented the analysis of pulsed cascaded paramet¬ 
ric down-conversion based on a monolithically integrated 
model system. For our real-world scenario, we derived 
the expectable two- and three-fold coincidence probabil¬ 
ities for the photon triplet detection process, where we 
included transmission and detection inefficiencies as well 
as noise contributions. 

Taking reasonable experimental apparatus into ac¬ 
count, we predicted photon triplet generation rates of 


~ 4 per hour, when the mean photon number of the 
primary PDC stage were (to) = 0.25 at laser repetition 
rates of 10 MHz. The impact of uncorrelated noise con¬ 
tributions to the signal-to-noise ratio of three-fold coin¬ 
cidences has been shown, and we inferred that a careful 
analysis of higher-order photon contributions is as essen¬ 
tial as the estimation of higher-order PDC processes in 
the first stage. 

Significant improvements of the detectable photon 
triplet rate become feasible using highly efficient detec¬ 
tors with very low dark count rates with efficiencies of 
over 80%. Likewise, the influence of higher-order pho¬ 
ton contributions on the primary idler detection rates 
will decrease using. Additionally, transition-edge detec¬ 
tor systems operate in free-running mode with intrinsic 
photon number resolution. Thus, a convenient separation 
of genuine photon triplets from higher-order three-fold 
coincidences becomes an attractive option. 

Higher repetition rates of the pump laser and the 
simultaneous reduction of the mean photon numbers 
per pump pulse are suitable options to furthermore 
reduce the impact of higher order photon contributions. 
Technological improvements to the individual PDC 
conversion efficiencies are also feasible using reverse- 
proton-exchanged waveguides [25] , 
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Appendix: List of realistic conditions for our model 
sytem 

In this appendix, we give an overview on reasonable 
numbers for the properties of the model system, i. e., 
individual detection efficiencies, conversion efficiencies, 
transmittances, and detector noise characteristics for the 
calculations in the main text. 
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TABLE I. Summary of assumptions for our calculations 

Optical element 

characteristic 

formula symbol 

transmittance 

pump 

pulse width 

Tpulse 

4.4 • 10" 11 s 

pump 

repetition rate 

-Rrep 

10 7 Hz 

waveguide 

endface transmittance 

??ef(A ~ 1590 nm) 

0.995 

idler 1 

noise probability per ns 

Pnoise,il 

7- 10 _t> 


detection efficiency 

»?Det,ii(A = 1625 nm) 

~ 0.45 


coupler efficiency 

^coupler,il (A — 1625 nm) 

~ 0.94 


waveguide transmittance 

WG.ii (A = 1625 nm) 

~ 0.92 


optical components transmittance 

J?oc,ii(A = 1625 nm) 

~ 0.3 


overall channel efficiency 

??ii 

~ 0.117 

signal 1 

coupler efficiency 

^coupler,si (A — T90.5 HUl) 

~ 0.995 


waveguide transmittance 

WG,si(A = 790.5 nm) 

~ 0.93 


2 nd PDC efficiency 

Dnom 

r PDC,2 

~ 2.7- 10 -7 


internal process efficiency 

PpDC,2 

1 

O 

CM 

lO 

CM 

signal 2/idler 2 

integral transmittance of optical components 

»?OC,s2/f?OC,i2 

~ 0.411/ ~ 0.292 


integral detector efficiency 

»?Det,s2/??Det,i2 

~ 0.25/ ~ 0.65 


noise probability per ns 

pnoise ,s2 //Aioise. i2 

7.5- 10 _6 /1.8 • 10 -5 


overall channel efficiency 

Vl2/V i2°' 

~ 0.103/ ~ 0.190 



